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A New Phenomenon: Sub-T^, Solid-State, 
Plasticity-Induced Bonding in Polymers 

Nikhil Padhye, David M. Parks, Bernhardt L. Trout, Alexander H. Sloeum 

Polymer self-adhesion due to the interdiffusion of macromolecules has heen 
an active area of research for several decades |[70l |4^ HU |72l |7^ |4l]|. Here, 
we report a new phenomenon of suh-Tg, solid-state, plasticity-induced bonding; 
where amorphous polymeric films were bonded together in a period of time on 
the order of a second in the soUd-state at ambient temperatures nearly 60 K below 
their glass transition temperature (Tg) by subjecting them to active plastic defor¬ 
mation. Despite the glassy regime, the bulk plastic deformation triggered the req¬ 
uisite molecular mobility of the polymer chains, causing interpenetration across 
the interfaces held in contact. Quantitative levels of adhesion and the morpholo¬ 
gies of the fractured interfaces vaUdated the sub-Tg, plasticity-induced, molecular 
mobilization causing bonding. No-bonding outcomes (i) during the compression 
of films in a near hydrostatic setting (which inhibited plastic flow) and (ii) between 
an ‘elastic’ and a ‘plastic’ film further established the explicit role of plastic de¬ 
formation in this newly reported sub-Tg solid-state bonding. 

If two pieces of a glassy polymer are brought into molecular proximity at temper¬ 
atures well below their glass transition temperature (Tg), negligible adhesion due to 
interdiffusion of macromolecules will be noted. Because polymer chains are kineti- 

*The supplementary videos can be obtained by emailing npdhye@mit.edu or npdhye@gmail.com, or can 
be viewed at http://web.mit.edu/npdhye/www/supplementary-videos.html 
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cally trapped well below the Tg E [301 EH 167), the time scales for relaxations in the 
glassy state are extremely large li22l[3^[J7l . Therefore, the system is completely frozen 
with respect to any cooperative segmental motions (a-like relaxation) |j71 that would 
cause interdiffusion. For example, the glass transition state itself is typically charac¬ 
terized by viscosity and diffusivity values of 10*^ Poise and 10^^^ m^/s, respectively, 
i6)- In El, assuming a viscosity of 10*^ Poise at the glass transition temperature, 
self-diffusion coefficients of forty polymers were estimated to be approximately 10^^^ 
m^/s (see Supplementary Section 3.1 for discussion). 

However, if the two pieces are brought into contact at a temperature above the 
glass transition temperature with the application of moderate contact pressure IjTQlEIl 
polymer chains from the two sides interdiffuse on experimental 
timescales. As a result of this interdiffusion, there is an optical disappearance of cracks 
and the development of strong bonds between the two surfaces over time. The strength 
of the developing interface is a function of temperature, time of healing and pressure, 
and the healing process continues until the interface acquires the bulk properties. Typ¬ 
ically, for times smaller than the reptation time, the interface toughness (Gc) and shear 
strength (o*) show a monotonic time-dependent growth as Ge ^ t^/^ and Oi ~ t'/^ 
li48l l2n l72l ISOi . The temperature strongly dictates the molecular mobility, with the 
self-diffusion coefficient of polymer melts usually ranging between and 10^^^ 

m^/s (see Supplementary Table 3). Moderate contact pressures (ranging from 0.1 MPa 
to 0.8 MPa) have been reported to be essential for facilitating the intimate contact be¬ 
tween the interfaces that allows interdiffusion. The chemical structure, the molecular 
weight and polydispersity of the polymer, the geometry of the joint, and the method of 
testing are critical factors affecting the measured strength of the interface. 

In the past two decades, there have been reports of polymer adhesion due to in¬ 
terdiffusion at temperatures below the bulk T^ with relatively long healing times of 
several minutes IMKIl, hours IE] da, and even up to a day ini. Such studies have 
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claimed that bonding due to interdiffusion at temperatures below the bulk Tg is possi¬ 
ble because the surface layer of a glassy polymer is in a rubbery state. The presence 
of a rubbery-like layer at the free surface with enhanced dynamics has been verified 
with experiments ll58l 132] and computer simulations ESI. The mean configuration of 
the macromolecules at the free surface is also perturbed in the direction normal to the 
surface. The resultant effect of entropic and enthalpic factors can lead to segregation 
or repulsion of chain ends at the free surface Ho). The segregation of chain ends at the 
free surface is also responsible for causing the depression of the glass transition tem¬ 
perature at the surface ioiisiiiisiiisi. However, such effects decay within distances 
comparable to the radius of gyration of the polymer. 

Although the motion of macromolecules in a glassy state is effectively frozen, 
stress-induced molecular mobility of glasses has been studied since the work of Eyring 
ED. Argon and co-workers demonstrated that the case II sorption rates of low 
molecular weight diluent species into a plastically deforming glassy poly(ether-imide) 
were dramatically enhanced, and were comparable with the sorption rates into the poly¬ 
mer at Tg, and that plastically deforming glassy polymers exhibit a mechanically di¬ 
lated state, which is representative of the molecular-level conformational rearrange¬ 
ments, such as those at Tg. A related study Q also reported an increase in the case II 
front velocity (of approximately 6.5 times) when an out-of-surface tensile stress was 
applied. Lee et al. ED showed that uniaxial deformation of PMMA 19 K below its 
Tg exhibited an increased molecular mobility by up to 1000 times. In E4l . the au¬ 
thors used NMR to probe deuterated semi-crystalline Nylon 6 and reported enhanced 
conformational dynamics in the amorphous regions of Nylon when deformation was 
carried out near Tg. Molecular dynamics simulations EOll also revealed increased tor¬ 
sional transition rates and thus enhanced molecular mobility during active deformation 
of a glass. The plastic deformation of glassy polymers is understood in terms of local¬ 
ized step-like shear cooperative displacements of lengthy chain segments, and the unit 
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plastic rearrangements are known as shear transformations il. According to molecu¬ 
lar dynamics simulations ll68ll . slippage of chains is the underlying feature of a shear 
transformation (for a detailed discussion, see Supplementary Section 3.2). Here, we 
report that active plastic deformation of glassy polymeric hlms held in intimate contact 
triggers requisite molecular-level rearrangement to cause interpenetration of polymer 
chains across the interface, which leads to bonding. Figure [T] shows the comparison 
of polymer self-adhesion through interdiffusion and the plasticity-induced technique 
proposed herein. 


'Interdiffusion' above Tg 



Figure 1: Comparison of polymer self-adhesion through diffusion at temperatures 
near or above Tg and the newly proposed suh-Tg, solid-state, plasticity-induced 
bonding in which bulk plastic deformation triggers the requisite molecular mobil¬ 
ity for chain interpenetration across the interfaces. 


Polymeric films were prepared by solvent casting (as shown in Figure using a 
base polymer (hydroxypropyl methylcellulose) and a plasticizer (polyethylene glycol, 
PEG-400). The base polymer HPMC was available under the trade name METHO- 
CEL in E3 and E15 grades. The molecular structures are shown in Eigure The films 
were assigned a name depending on the base polymer and weight percent (wt.%) of the 
plasticizer in the film with respect to the base polymer (see Methods and Supplemen- 
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Polymer solution 


Solvent casting 


Polymer film after drying 


Figure 2: Steps involved in the preparation of polymer films through solvent cast¬ 
ing: (1) a homogeneous solution of polymer and plasticizer in ethanol and water, 
(2) spreading of the solution on a glass surface via a knife, and (3) evaporation of 
solvents and formation of a glassy film after drying. 





HPMC 


PEG 


Figure 3: Molecular structures of hydroxypropyl methylcellulose (HPMC) and 
polyethylene glycol (PEG). 


C \ 



Figure 4: True stress-strain curves for three film formulations: E3-alone-42.3% 
PEG, E3/E15 in 1:1-42.3% PEG and E15-alone-42.3% PEG at ambient tempera¬ 
tures. The nominal strain rate for tensile testing was chosen as 0.0025 sec '. 













tary Section 1). Films made from E3-alone-42.3% PEG, E3/E15 in 1:1-42.3% PEG 
and E15-alone-42.3% PEG exhibited values in the range of 72-78°C. (See Methods 
and Supplementary Section 2.2). Their true stress-strain curves in tension are shown 
in Eigure All three hlms exhibited ductility, represented by their ability to undergo 
plastic flow. 



Figure 5: Roll-bonding was achieved by passing a stack of film layers with a total 
initial thickness ti between compression rollers to yield a final-thickness t 2 . The 
peel-test was carried out on roll-bonded sample at the middle interface. 



Figure 6: Lap specimens were prepared between two film layers by applying 
compression loads on the overlapping area. Lap shear-strength measurements 
were performed in a tensile mode. 


Bonding experiments were carried out at ambient conditions, (i) Stacks of six film 
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layers (each layer ^100 /rm) were fed through a roll-bonding machine to achieve active 
plastic deformation at ambient temperatures. Peel tests were performed to measure the 
mode I fracture toughness (Ge [J/m^ ]), Figure and (ii) lap specimens were prepared 
to measure the shear-strength (a^ [MPa]), Figure |^(see Methods and Supplementary 
Section 4 for details on roll-bonding, peel testing and lap shear strength testing). Gc 
represents the work done per unit area for debonding the interface during a peel test. 
Ci indicates the maximum shear stress sustained by the bonded interface before failure. 
The effective thickness reduction was used as a measure of plastic strain during bonding 
in all of the cases. 

Figure j^shows a snapshot of several layers of the him (E3/E15 in 1:1-42.3% PEG) 
with an initial thickness of ti=0.60 mm undergoing roll bonding through active plastic 
deformation with a hnal thickness reduced to I 2 =0.533 mm (see supplementary video 
SI). 



Stack of film layers Multiple layers undergoing sub-Tg, solid-state, roll-bonding Roll-bonded layers 

-litial thickness 0.024" {0.60 mm) (final thickness 0.021" (0.53 mm) 


Figure 7: Illustration of sub-Tg, solid-state, plasticity-induced roll bonding of 
E3/E15 in 1:1-42.3% PEG films nearly 60 K below Tg. For this case, the nominal 
thickness strain is ep= |f 2 — 1 \ |/fi=11.7%. 

The Ge results for the three films are shown in Figure G^ correlated with the 
plastic strain in a non-monotonic fashion, first increasing and then decreasing. The ad¬ 
hesion between two interfaces held together by van der Waals forces, hydrogen bonds, 
or chemical bonds can only give Gc values in the range of 0.05 J/m^ , 0.1 J/m^ and 
1.0 J/m^, respectively ll27ll . The surface energy of glassy polymers itself is quite small 
lEl (on the order of 0.08 J/m^); therefore, negligible adhesion is noted when two such 
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Figure 8: Fracture toughness (Ge [J/m^]) versus plastic strain plots for E3/E15 in 
1:1-42.3%PEG, E3-alone-42.3%PEG and E15-alone-42.3%PEG. 


0.14 

0.12 

0.1 

0.08 

0.06 

0.04 

0.02 


E3-alone-42.3% PEG 
E3/E15 in 1:1>42.3%PEG 
E15>alone-42.3% PEG 


A 

A 


A ^ O 
.O O 


O 


o 


o 


o 


o 


o 


8 10 12 
Plastic Strain (%) 


CA 


Figure 9: Lap shear-strength (c^ [MPa]) versus plastic strain plots for E3/E15 in 
1:1-42.3%PEG, E3-alone-42.3%PEG and E15-alone-42.3%PEG. 


surfaces are brought into mere molecular proximity. However, glassy polymers can 
exhibit higher fracture toughness owing to the irreversible deformation of the macro¬ 
molecules. The quantitative levels of Gc obtained here, with a maximum value nearly 
10 J/m^, could only be attributed to the irreversible processes of chain pull-outs, disen¬ 
tanglement and/or scissions during debonding, which could only happen if plasticity- 
induced molecular mobilization and chain-interpenetration led to bonding. Even poly¬ 
mer adhesion leading to Gc values as low as 1.2 J/m^ ifT^ and 2.0 J/m^ IItTII has been 
attributed to irreversible chain pull-out mechanisms during fracture. Other mechanisms 
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of adhesion such as acid-base interactions, capillary effects, electrostatic forces and/or 
any other conceivable mechanism do not apply in the current context (for a detailed dis¬ 
cussion on the types of forces giving rise to adhesion, see 153]). The shear strength (Oi) 
plots, shown in Figure also exhibited a non-monotonic correlation with the bonding 
plastic strain. Quantitative levels of the values reported here compare with those 
in m, in which adhesion due to interdiffusion of chains below the bulk over long 
times, on the order of several minutes, was reported. The reported levels of bulk plastic 
strains also rule out mechanical interlocking of asperities to cause adhesion because, 
at the levels of plastic strains reported here, the surface asperities would necessarily 
flatten out. Surface characterization of the hlms through AFM, before bonding, re¬ 
vealed nano-scale roughness (Ra) on the order of 6.91-22.7 nm (see Supplementary 
Section 2.6). By contrast, high levels of plastic strain led to an increased contact area, 
and if factors other than chain interpenetration were responsible for bonding, we would 
expect a monotonic increase in Ge or Cj. The lowering of Gc or Os at high levels of 
plastic strain could be explained on the basis of anisotropic growth in the microstruc¬ 
ture such that the polymer chains oriented in the direction of the principal stretches 
(compression and rolling directions). We suggest that increasing levels of such chain 
orientation ultimately lead to less effective chain interpenetration across the interface, 
which diminishes bonding at higher strain. 

A comparison of the surface morphology before bonding and after the fracture is 


shown in Figure 10 The debonded surfaces indicated events of chain scissions or pull¬ 
outs due to fracture, which were similar to those reported upon fracture of polymers 
welded through interdiffusion 117211^ [141 . 

To explicitly demonstrate the role of bulk plastic deformation, we designed a ‘hy¬ 
drostatic die’ setup, which was capable of generating high levels of hydrostatic pressure 


while inhibiting the macroscopic plastic flow. Figure 11 shows a comparison in which 
a stack of hlms (E3/E15 in 1:1-42.3% PEG) was compressed (i) without any constraints 
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Fractured surfaces after de-bonding 


Figure 10: SEM images of E3/E15 in 1:1-42.3% PEG, E3-alone-42.3% PEG, and 
E15-alone-42.3% PEG, films before bonding and after debonding. The nominal 
plastic strains during roll-bonding for E3/E15 in 1:1-42.3% PEG, E3-alone-42.3% 
PEG, and E15-alone-42.3% PEG, films were 15.53%, 8.12%, and 10.18%, respec¬ 
tively. 


and (ii) with the ‘hydrostatic die’ constraint. In the first case, the stack underwent 
macroscopic plastic flow, and the layers bonded to form an integral structure, whereas 
in the case of the hydrostatic die constraint, no permanent thickness change was ob¬ 
served, and the layers easily splayed apart after removal (see Supplementary video S2). 
In another experiment, we attempted to roll-bond E3/E15 in a 1:1-0% PEG film with 
E3/E15 in a 1:1-42.3% PEG film. Eilms with 0% PEG exhibited negligible plastic flow, 
and were therefore incapable of being subjected to plasticity-induced molecular mobi¬ 
lization, which also led to a no-bonding result (see Supplementary video S3). Nanoin¬ 
dentation experiments clearly revealed the differences between the elastic response of 
the 0% PEG films and the plasticity of the 42.3% PEG films (see Supplementary Eig- 
ure 16). Both experiments demonstrated that activating bulk plastic flow on both sides 
of the interface was an essential requirement for bonding. This result also strongly 
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stack of film layers Compression of layers: "Simple upsetting" Bonded layers 

(initial thickness 0.84 mm) (final thickness 0.70 mm) 


(a) 



Stack of film layers inside die 

(initial thickness 0.84 mm) 


Compression inside die 


After compression 

(no measurable change in thickness) 


Layers Splay 


(b) 

Figure 11: Compression of stacks of films (a) without any die containment to per¬ 
mit macroscopic plastic How and bonding, (b) in a hydrostatic die that is capable 
of generating high levels of hydrostatic pressure but limits the plastic fiow, and 
consequently no bonding takes place. In both cases peak nominal compressive 
stresses (78.98 MPa) was kept same. 


reflected that effects such as the presence of a rubber-like layer at the surface, where 
the Tg may be lower than the bulk Tg and the time scales for segmental relaxations 
may be relatively small, by itself could not lead to adhesion of the magnitude observed 
during roll-contacts lasting on the order of a second (see supplementary Table 4 for 
estimates of the rolling times). Bonding below the bulk Tg (without any bulk plastic 
deformation), as reported in the literature, requires substantially longer durations. Fur¬ 
thermore, the existence of any enhanced relaxation of the polymer chains (or segments) 
in the surface layer would be severely restricted by any portions of the macromolecules 
extending into the glassy-bulk beneath; hence, long-range diffusion in a short time is 
not possible. Finally, although not considered in these prior reports, it is plausible that 
moderate contact pressures, applied over relatively long healing times at temperatures 
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near T^, contributed to mechanically enhanced molecular mobility that led to bonding 
via mechanisms similar to those described here. 

Although rapid plastic deformation can cause a temperature increase, fully adia¬ 
batic analysis revealed an upper bound temperature increase of only 3.6°C (see Supple¬ 
mentary Section 3.3). The mechanically activated polymer mobility is mechanistically 
quite different from molecular mobility at temperatures above Tg. The shear transfor¬ 
mation units of plastic deformation are also accompanied by local transient dilatations 
(volume changes) that can facilitate opportunities for establishing entanglements across 
the interface so that bonding can take place on the order of a second. The self-diffusion 
coefficient (D) of a polymer chain in its melt state shows a strong dependence on the 
molecular weight, D ^ M or D ~ in accordance with the Rouse or the reptation 
model, respectively. However, all three blends of polymer considered here, E3-alone, 
E15-alone and E3/E15 in 1:1, were roll-bonded on the order of a second, which was 
a significant contrast from the mechanism of polymer adhesion due to interdiffusion. 
The unprecedented aspects of the newly reported phenomenon and underlying mecha¬ 
nisms are expected to open new avenues for research and applications. 

Methods 

Film-Making: 

Hydroxypropyl methyl cellulose (HPMC), trade name METHOCEL, in grades E3 and E15 was obtained 
from Dow Chemical (Midland, Michigan, North America). PEG-400 was purchased from Sigma-Aldrich 
(Milwaukee, Wisconsin, North America). Appropriate amounts of E3, E15 and PEG were mixed in desired 
amounts with ethanol and water, and a homogeneous solution was obtained through mixing with an electric 
stirrer for 24 h. After completion of the blending process, the solution was carefully stored in glass bottles at 
rest for 12 h to eliminate air bubbles. Solvent casting was carried out using a casting knife applicator from 
Elcometer (Rochester Hills, Michigan, North America) on heat-resistant borosilicate glass. All of the steps 
were carried out in a chemical laboratory where ambient conditions of 18° ± 2°C and R.H. 20%±5% were 
noted. The residual moisture content in the films after drying was measured using Karl Fischer titration. 

Bonding Experiments: 
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Roll bonding was carried out on a machine capable of exerting the desired load levels to achieve active 
plastic deformation. The 200 mm diameter rollers were operated at an angular speed of 0.5 rev/min, leading 
to an exit speed of 5.23 mm/s. Peel tests were caiTied out to measure mode I fracture toughness (see Sup¬ 
plementary video S4). Lap specimens were prepared using compression platens on an Instron mechanical 
tester. For both roll-bonded and lap specimens, for the sake of consistency, the adhesion measurements were 
canied out on the bonded interfaces between the top-top surfaces (exposed side during diying). Film layers 
were stacked accordingly. Top-bottom and bottom-bottom joining led to similar bonding results. The hydro¬ 
static die and upsetting experiments were carried out on the Instron. The roll-bonding machine and fixture 
for the peel test were designed and fabricated in Massachusetts Institute of Technology (Cambridge, North 
America) (see Supplementary Section 4). 

Characterization: 

The moleculai* weights of E3 and E15 were estimated from viscosity measurements. The amorphous 
nature of the films were verified by XRD. SEM and AFM were performed to analyze the surfaces. DMA 
was performed to determine the Tg. Tensile stress-strain curve tests, fracture toughness through peel tests, 
and lap shear tests were carried out. Nanoindentation was carried out to measure the hardness. The specific 
heat capacity was measured using DSC. 

The X-ray diffraction patterns were recorded using a PANalytical X’Pert PRO Theta/Theta powder X- 
ray diffraction system with a Cu tube and an X’Celerator high-speed detector. AFM images were obtained 
using a Dimension 3100 XY closed loop scanner (Nanoscope IV, VEECO) equipped with NanoMan soft¬ 
ware. Height and phase images were obtained in tapping mode in ambient air with silicon tips (VEECO). 
DMA was carried out on a TA Q800 instrument. Mechanical testing was performed on an Instron mechan¬ 
ical tester. Nanoindentation tests were canied out on a Triboindenter Hysitron instrument. Calorimetry was 
performed on a TA Q200 instrument. The viscosity was measured on an HR-3 Hybrid rheometer. 

Acknowledgements: 

The authors acknowledge the funding from Novartis Pharma AG and facilities at the Novartis-MIT 
Center for Continuous Manufacturing program ED where this research was carried out. N.P. also appreciates 
discussions with Professor Robert E. Cohen. 

Author Contributions: 

N.P. conducted experiments, designed and developed the experimental set-ups, and wrote the letter 
and supplementary. D.M.P. identified the sub-T^, solid-state, plasticity-induced bonding and provided su¬ 
pervision on the subject of mechanics. B.L.T. advised on the phenomenological aspects of adhesion and 
supervised the overall development of thin-film technology at MIT. A.H.S. supervised the design and devel¬ 
opment of mechanical fixtures and machines. D.M.P, B.L.T. and A.H.S jointly supervised the work. All 


13 


authors contributed to review of the manuscript, and participated in discussions during this research. 


14 



1 SUPPLEMENTARY FILM-MAKING 


Table 1; Formulations employed in making polymer films from HPMC E3 and E15 
with different levels of plasticizer (the amounts have been rounded off to nearest 
grams). 


Polymer film 

Composition 


E3 

E15 

Water 

EtOH 

PEG 


(g) 

(g) 

(g) 

(g) 

(g) 

E3/E15in 1:1-0% PEG 

15 

15 

96 

96 

0 

E3/E15in 1:1-28.5% PEG 

15 

15 

96 

96 

12 

E3/E15 in 1:1-42.3% PEG 

15 

15 

96 

96 

22 

E3/E15in 1:1-59.5% PEG 

15 

15 

96 

96 

44 

E3-alone-42.3% PEG 

30 

0 

96 

96 

22 

E15-alone-42.3% PEG 

0 

30 

96 

96 

22 


Polymeric films employed in this study comprise of a base polymer hydroxypropyl 
methyl cellulose (HPMC) and a compatible plasticizer, PEG-400. HPMC is cellulose 
ether and available from Dow Chemical under the trade name of METHOCEL. We 
acquired METHOCEL E3 and METHOCEL E15 products from Dow Chemical. 

HPMC is an uncrosslinked polymer and shows an excellent film formability due to 
its underlying cellulose structure in which all the functional groups lie in the equatorial 
positions, causing the molecular chain of cellulose to extend in a more-or-less straight 
line and easing the formation of the film. Table shows the sample weights of the 
contents used in preparation of the solutions. Eilms were made through casting and 
drying of the prepared solutions. As seen in Table [T] films are referred to based on the 
amounts of E3, E15 and Wt.% of PEG-400. Eor example, E3/E15 in 1:1-42.3% PEG 
implies that E3 and E15 are present in one-to-one ratio and the Wt.% of PEG in the 
film is 42.3%, since 22 g PEG in 15 g E3 plus 15 g E15 is 42.3%. 

Karl Eischer titration was carried out to determine the residual moisture content in 
the films after drying. Dimethyl sulfoxide (DMSO) was used as a reagent for dissolving 
films. DMSO was purchased from Sigma Aldrich (ACS reagent grade). The solution 
of the dissolved film in DMSO was fed into the Karl Eischer Titrator, and the residual 
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Table 2: Residual moisture in films after drying, measured through Karl Fischer 
titration. %Wt. indicates residual moisture in the films after drying. 


Polymer film 

Residual H 2 O 

(%Wt.) 

E3/E15 in 1:1-0% PEG 

3.70 

E3/E15in 1:1-28.5% PEG 

7.21 

E3/E15in 1:1-42.3% PEG 

4.29 

E3/E15in 1:1-59.5% PEG 

2.45 

E3-alone-42.3% PEG 

2.92 

E15-alone-42.3% PEG 

4.54 


moisture was estimated. Table shows the average amounts (repeated three times) of 
the estimated residual moisture contents in the films. 


2 SUPPLEMENTARY CHARACTERIZATION 


2.1 Role of Plasticizer and Mechanical Properties 

PEG-400 acts as a compatible plasticizer for HPMC. Inclusion of a plasticizer within 
the polymer matrix enhances its ductility (or plastic flow characteristics). Plasticizers 
work by embedding themselves between the chains of polymers and spacing them apart 
by increasing the free volume. By dissolving and mixing intimately, PEG molecules 


disrupt the secondary bonds between the polymer chains. Eigure 12(a) illustrates this 
effect. 

Eigure |12(b)| shows true stress-strain behavior of films made from E3/E15 in 1;1 
with varying the PEG concentrations 0%, 28.5%, 42.3% and 59.5%, respectively. The 
plasticization effect of increasing PEG (Wt.%) is evidenced by the lowering of the 
initial modulus and the yield strength and, increase in the failure strain. The maximum 
failure strain occurs for 42.3% PEG film. Clearly, all films containing PEG demonstrate 
large ductility that is absent in 0% PEG film. The effect of including PEG on the 
lowering of glass transition temperature is discussed next. 
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strain 


(b) 

Figure 12: (a) Schematic role of plasticizer on molecular configurations, (b) Ef¬ 
fect of PEG-400 on tensile true stress-strain behavior of polymeric films. The 
tensile tests were carried out at ambient temperature and a nominal strain rate of 
0.0025 s 

2.2 Dynamic Mechanical Analysis 

Dynamic Mechanical Analysis was performed on all the films listed in Table [T] A 
temperature sweep was performed at 1 Hz frequency. Figures |13(a)| to |13(f)| show 
the plots of loss modulus, storage modulus and tan 5 for six different films. The glass 
transition temperature is determined by the peak in the tan 5. For E3-alone-42.3% PEG 
Tg was estimated to be 72° C, and for E15-alone-42.3% PEG and E3/E15 in 1:1-42.3% 
PEG Tg was estimated as 78° C. Inclusion of PEG evidently lowers the glass transition 
temperature and broadens the temperature range over which the glass transition takes 
place. 
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(a) E3/E15 in 1:1-0% PEG (b) E3/E15 in 1:1-28.5% PEG 






Figure 13: Dynamic Mechanical Analysis Curves. 

2.3 Molecular Weight 

Viscosity measurements for 2% aqueous solution of E3 and E15 were carried out based 
on the procedure prescribed by Dow El. The viscosity curves for E3 and E15 are 
shown in the Eigure E3 solution shows negligible rate dependence; with the vis¬ 
cosity lying in the range 3-4 mPa-s. E15 solution shows some rate dependence with 
viscosity lying in the range 14-21 mPa-s. Erom the Dow Chemical manual El. the 
viscosity for 2% aqueous solution of different grades E3, E5, E6, E15 and E50 are 
specified as 3, 5, 6, 15 and 50 mPa-s, respectively. In another Dow manual IT], the 
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Figure 14: Viscosity curves of 2% aqueous solution of METHOCEL-E3 and 
METHOCEL-E15. 


ranges for the viscosity of 2% Wt. solution of E3 and E15 are specified as 2.4 - 3.6 
mPa-s and 12-18 mPa-s, respectively. Our measurements overlap well within these 
specifications. If we choose a representative viscosity of 3.8 mPa-s for E3 and 16 mPa- 
s for E15, then based on the viscosity and molecular weight relationship from 121, we 
estimate the number average molecular weight (M„) for E3 and El5 approximately as 
8,200 and 20,000, respectively. 

A detailed molecular characterization of METHOCEL cellulose ethers presented in 
ISl, also led to estimation of weight average (M„,) and number average (M„) molecular 
weights as: (i) E3: M„ = 8,100 and = 20,300 with Mn./M„ = 2.5, and (ii) E15: 
M„ = 24,800 and = 60,300 with Mvv/M„ = 2.4. Such estimations are consistent 
with those we obtained. In the same study ll44ll . the degree of polymerization (DP) was 
reported as: (i) E3, DP= 77, and (ii) E15, DP=296, and the weight average radius of 
gyration (Rgw) as: (i) E3, Rg* = 7.4 nm, and (ii) E15, Rg* = 15.1 nm. 


2.4 X-Ray Diffraction 


HPMC is a cellulose derivative and well known to exist in amorphous form. To verify 
the amorphous characteristic, an X-ray diffraction was carried out on E3/E15 in 1:1- 
42.3% PEG, shown in Eigure As expected, a diffused pattern without any peaks is 
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Figure 15: XRD of E3/E15 in 1:1-42.3% PEG. 

obtained, thus, indicating absence of any crystallinity. 


2.5 Nanoindentation 



Figure 16: Nanoidentation load versus displacement curves for E3/E15 in 1:1- 
0% PEG, E3/E15 in 1:1-42.3% PEG, E3-alone-42.3% PEG and E15-alone-42.3% 
PEG films. Indentation experiments were carried out in load controlled mode 
with chosen peak loads up to 300 and 2000 for films with 42.3% PEG and 
0% PEG, respectively. 
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Nanoindentation experiments were performed on E3/E15 in 1:1-0% PEG, E3/E15 
in 1:1-42.3% PEG, E3-alone-42.3% PEG and E15-alone-42.3% PEG films. The ex¬ 
periments were carried out in a force controlled mode with a maximum force of 2000 
pN and 300 pN for 0% PEG and 42.3% PEG films, respectively. A larger load for 
the 0% PEG film was chosen in order to activate sufficient plastic indentation so that 
its hardness could be measured. Berkovich indenter with a root radius of 150 nm was 
used. The load versus displacement curves for all the films are shown in the Eigure 


16 The film with 0% PEG shows a relatively large indentation force and large elastic 


recovery, whereas films with 42.3% PEG films show little elastic recovery and large 
residual indentation depth. Based on these behaviors, the 0% PEG film can be called 
an ‘elastic’ film and the 42.3% PEG film as a ‘plastic’ film. 

Using Oliver-Pharr method we estimated the hardness from the nano-indentation 
tests. The hardness values for E3/E15 in 1:1-0% PEG, E3/E15 in 1:1-42.3% PEG, 
E3-alone-42.3% PEG and E15-alone-42.3% PEG films were 144.0 MPa, 10.83 MPa, 
10.151 MPa, and 11.48 MPa, respectively. This shows that the film with 0% PEG 
is “hard” and unlikely to demonstrate plasticity-induced molecular mobilization for 
bonding at the load levels where 42.3% PEG films exhibit sufficient plastic flow. This 
is also consistent with the no-bonding outcome between 0% PEG film and 42.3% PEG 
film, as shown in the Supplementary video S3. 


2.6 Atomic Force Microscopy 

Eigure shows sample AEM scans of a 5 pm x 5 pm area on the top surface of three 

L" lv-| 

films with 42.3% PEG, along with the average roughness given by Ra = y . The top 
surfaces of films exhibit nano-scale roughness, however, this scale of roughness does 
not play any important role when we have reported bulk plastic strains essential for 
bonding. 
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Height 1.0 pm Height 1.0 pm Height 1.0 pm 


(a) Top surface of E3/E15 in (b) Top surface of E3-alone- (c) Top surface of E15-alone- 
1:1-42.3% PEG film, Ra - 6.91 42.3% PEG film, Ra - 22.7 nm 42.3% PEG film, Ra - 8.63 nm 
nm 

Figure 17: Measurement of Nanoroughness using Atomic Force Microscopy. 

3 SUPPLEMENTARY DISCUSSION 

3.1 Polymer Dynamics and Self-Diffusion 

Polymer melts are an equilibrium system and their mobility is commonly described by 
the models of Rouse, reptation, etc. The fundamental dynamic property that character¬ 
izes the average motion of a polymer chain is the coefficient of self-diffusion (D). In 
the Rouse regime, polymer chains are unentangled and the diffusion coefficient Drouse 
where N is number of monomers. In the reptation regime polymer chains ex¬ 
perience topological constraints, and the lateral displacements in the entangled chain 
stracture are negligible compared with the longitudinal diffusion or “reptation” in the 
“tube” formed by neighboring molecules with Dreptation The transition from 

Rouse to reptation regime is often noted at the critical entanglement length Ng. The 
self-diffusion coefficient in polymer melts strongly depends on the temperature, molec¬ 
ular weight, and molecular characteristics. Typical values of self-diffusion coefficient 
are listed in Table 

Usually, the diffusivities of the polymer melts are quite small compared to the liq¬ 
uids composed of simple molecules (for example D„ater at 298.1 K and 1 atm is nearly 
2.3 X10^® m^/s 1491 ). However, high diffusion coefficients for polymer melts can be 
noted when molecular weights are small and temperatures are well above the melting 
point or the glass transition temperature. 
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Table 3: A Short Summary of Self-Diffusion Coefficient (D) from literature. T„, 
and stand for melting and glass transition temperature, respectively. 


Ref. 

Polymer 

Mol. Wt. 

(g/mol) 

Temp. 

(K) 

D 

(m^/s) 

Comments 

GIl 

Linear H-PB 

5x10'*-20 X 10'* 

398.15 

10-*'*-5x10-"' 

T„, 

~ 381.15 K 

(Hill 

Polyisoprene 

560-9.82X 10'* 

373.15 

2.2x10-*°-1.0x10-*'* 

T* 

~ 192X0 

(H 

Polybutadiene 

69Q-4.99xW‘^ 

373.15 

7.0x10-** -2.0x10-*'* 

Ts 

< 183.15 K 

(SD 

Polyethylene 

200-12x10'* 

448.15 

6.6x10-*°-1.3x10-*'* 

T„, 

-353.15-400.15 K 

ED 

PDMS 

500-5x10^ 

293.65 

7.0x10-*°-5x 10-*^ 

Ts 

- 150.15 K 

(H 

PS 

600- 19x10^ 

487.8 

2.5x10-*°-1.5x 10-*° 


- 333.15-373.15 K 


When the temperature of a glass forming liquid is lowered and the glass transition 
temperature is approached from the above, kinetics of a glass-forming system shows a 
drastic slow down, and time scales for relaxations become orders of magnitude larger 
than at higher temperatures representative of the melt. As an example, in Il69l the self- 
diffusivity of a single component glass former, tris- naphthylbenzene (TNB), falls from 
D = 10^'^ m^/s at 405 K to D = as Tg (close to 335 K) is approached. In ll5^ it 
was shown that as Tg is approached from above the diffusivity for o -terphenyl (OTP) 
is of the order m^/sec. Similar slow downs of self-diffusion, with D approach¬ 
ing m^/s near the glass transition temperature, were noted for o-terphenyl and 

metallic melts of PdCuNiP systems 16^ . 

The Bueche-Cashin-Debye equation ifTOlfTSll . which relates diffusivity and viscos¬ 
ity, is given as: 


Dri AKbT 

p 36 M ^ 

In the above equation, A is the Avogadro constant, Kb is the Boltzmann constant, 
T is the absolute temperature, is the mean-square end-to-end distance of a single 
polymer chain, and M is the molecular weight. If we estimate D for our polymer at the 
glass-transition temperature by considering r|=10'^ Poise, p=l 180 Kg/m^, R^= 6x7.4^ 
nm^ (using Rg of E3, and R^=6x Rg) and M=20,300 g/mol, T = 352 K, the estimated 
value of D is 1.12x 10^^^ m^/s. This is a remarkable estimate in terms of the order 
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of magnitude and compares well with the self-diffusivity of 10^^^ m^/s as reported in 
II 52 I . Such small diffusivity near Tg clearly indicates the aspect of kinetic arrest as the 
glass transition temperature is approached. 

If we consider a scenario; in which a diffusion distance of x=10 nm is to be 
achieved in a time of one second, then D must be greater than 0.5 X 10 m^/s. This 

is impossible in the solid-state, 60 K below the bulk-Tg, and clarihes the distinction of 
polymer welding above Tg with respect to newly reported plasticity-induced molecular 
mobilization and bonding which occurs in a period of time on the order of a second. 


3.2 Stress-Induced Molecular Mobility and Plastic-Deformation 

In 1936, Eyring undertook an absolute reaction-rate approach to demonstrate the effect 
of stress on lowering of viscosity in certain gels, glasses or crystals and enhancement 
of molecular mobility. According to Erying’s model; 


1 

Vl 


Ta(7’,T) =lo-exp 


Ea-tV* 

KbT 


Here, Vi is the jump frequency (inverse of the relaxation time) of the molecules in 
the direction of applied stress. To is a constant (vibration time scale), T is the applied 
shear stress, V* is the activated volume, Ea is the potential barrier that molecules have 
to overcome in going from one conhguration to another. The above equation states that 
the relaxation time decreases in the direction of the applied shear-stresses and molecu¬ 
lar transport is facilitated in the direction of shear stresses. However, Eyring’s model is 
found to be applicable only in the regimes of linear visco-elasticity for polymers, and 
fails to capture the dramatic changes in the energy landscapes associated with plastic 
deformation of a glass. 

As stated in the main letter, plastic deformation in polymers at a continuum scale 
is understood in terms of shear transformations i.e. events of spatial rearrangements of 
molecular clusters causing stress-relaxation. Consider the scenario shown in the Eigure 
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(a) 




Figure 18: Mechanism of plastic deformation and shear transformation in glassy 
polymers (a) A unit shear transformation in a kinetically trapped state under 
shear stress comprises of an initial elastic shear-strain which is followed hy plastic- 
relaxation of polymer chain segments, (h) Free energy landscape, associated with 
a polymer chain, during a unit shear transformation (c) Accumulation of several 
shear transformations leads to macroscopic plastic deformation. 
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well below the polymer chains are kinetically trapped in their local configura¬ 


tions and timescales for mobility (specifically translation motions) of these chains are 


extremely large. However, application of shear-stress on the material element causes 
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its deformation and the polymer chain under consideration changes its orientation; first 
elastically, and then due to some local perturbation it relaxes plastically while over¬ 
coming the potential barrier set up due to neighboring molecules. 

If no stresses were applied and temperature were held far below then the transi¬ 
tion of the mean configuration of a polymer chain from its kinetically trapped configu¬ 
ration would not happen on experimental time scales. However, qualitatively speaking, 
the application of stress enhances the mobility of the polymer chain as it relaxes, and 
changes its configuration on experimental time scales. Effect of such cumulative events 
during active plastic deformation characterize the enhanced dynamics in deforming 
glasses below Tg. 

We emphasize that by kinetically trapped state of a glass it is implied that any co¬ 
operative segmental relaxations or long range diffusive motions of chains are severely 
restricted; however, secondary relaxation processes (those corresponding to vibrations 
of side groups like p, y, 5, etc. relaxations) may still be active. But, such weak sec¬ 
ondary relaxation processes are incapable of giving any interdiffusion and pronounced 
adhesion in a short-time when two interfaces are brought together in molecular prox¬ 
imity, unless enhanced mobility is triggered through plastic deformation. 

The fundamental differences between polymer mobility at high temperatures (well 
above Tg) and stress-assisted molecular mobility (well below T^) can be summarized 
as follows: The motion of polymer chains (or segments) in a polymer melt can be 
described based on diffusion models. It primarily occurs due to high kinetic energy of 
the polymer chains (or segments), and available free-volume (or physical space) due to 
which chains (or segments) can sample new orientations effectively. The polymer melts 
(above T^) are spatially homogeneous and thermodynamically in an equilibrium state, 
whereas, plastic deformation and associated enhanced mobility in a glassy polymer 
is not at all an equilibrium concept. The root mean square displacement of center of 
mass of a polymer chain will increase monotonically with time during diffusion in a 


26 



polymer melt, however, the mechanically assisted enhanced mobility in polymers only 
occurs during active plastic deformation and stops when plastic straining stops. The 
average kinetic energy of a polymer molecule is large in a polymer melt compared 
to that in the solid-state glass well below T^. Typical values of activation energy for 
diffusion in molecular liquids at room temperature can be as low as 5-10 kcal/mol and 
therefore diffusion can be thermally activated. However, for a shear-transformation the 
activation energy (for example inorganic glasses is 350-400 kcal/mol [SI) and therefore 
plastic deformation is not thermally activated at room temperatures on experimental 
time scales. Although plastic deformation can be accompanied by a temperature rise, 
at relatively slow strain-rates the associated temperature rise is negligible. 


3.3 Temperature Rise 

It is worth checking for any temperature rise due to irreversible mechanical work during 
plastic deformation, and if such temperature rise is responsible for enhanced molecu¬ 
lar mobility leading to bonding. As an illustration, we measured the specific heat of 
E3/E15 in 1:1-42.3%PEG film through differential scanning calorimetry, as shown in 


Eigure 19 Erom the rate of heat flow into the sample and specified rate of temperature 
rise during thermal scan the Cp is obtained as 1860 J/Kg-K, and the density was mea¬ 
sured to be p = 1180 Kg/m^. Based on the stress-strain curves if we estimate the flow 
stress for plastic deformation to be = 8 MPa, then for a plastic strain of Ep = 0.5, 
the adiabatic temperature rise is estimated to be: 


AT = 3.6°C 

pCp 

As seen here, the temperature rise according to fully adiabatic analysis is quite 
small. External work due to the application of stresses leads to mechanically-assisted 
(and not temperature assisted) enhanced mobility of polymer chains (or segments). The 
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operative micro-mechanisms of plastic-relaxation at a molecular level are dependent 
on particular molecular characteristics, and may at best be explored through computer 
simulations. In our opinion, this remains as an open question. 
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Figure 19: DSC scan of E3/E15 in 1:1-42.3% PEG film. 


4 Bonding Experiments 


4.1 Roll-Bonding Machine 


Figure 20 shows a CAD model of the roll-bonding machine designed for this work. 
The machine is capable of achieving different levels of plastic strain by adjustment 
of the gap between the rollers and monitoring the compression load during rolling. 
The angular speed of the rollers is controlled using a stepper motor. The radius of the 
rollers (R) is 100 mm, much larger than the total initial thickness of a film-stack (ti), 
which is typically less than 1 mm. The incoming stack of film behaves like a thin strip 
and through-thickness plastic deformation is triggered under such conditions. From 
kinematics of rigid-plastic rolling of thin-strip lIMl the time spent during active plastic 
deformation can be estimated as follows: 


T = 


-ti) 

V2 


( 2 ) 
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Motor 


Figure 20: A roll-bonding machine to carry out sub-Tg, solid-state, plasticity- 
induced bonding. 

In the above equation, t] is the initial thickness of film-stack, t 2 is the thickness of 
film-stack at the exit, V 2 is the linear speed at the exit. For V 2 = 5.23 mm/s, ti=.6 mm 
and t2=0.45 mm (indicating 25% nominal plastic strain), the time spent by a material 
element under the roller would be approximately 0.74 s. This is how we achieve sub- 
Tg, solid-state, plasticity-induced roll-bonding in a period of time on the order of a 
second. The Supplementary video S1 demonstrates how a stack of films with a certain 
initial thickness is subjected to active plastic straining leading to sub-Tg, solid-state, 
plasticity-induced bonding. The final thickness of the roll-bonded stack is less than the 
initial. Complete details on the roll-bonding machine and process will be available in 

63. 


4.2 Mechanics of Peel Test 

Figure j^shows a snapshot of the peel test. A peel test fixture was designed to perform 
accurate mode-I fracture. Such a test is also commonly known as T-peel test in the lit- 
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Figure 21: Peel-Test in mechanical tester to determine mode-I fracture toughness 
(G,). 


/-\ 



Displacement (mm) 

V_ ! _/ 


Figure 22: Force versus displacement curve during Peel-Test. In the steady-state 
peeling the peel-force becomes steady with respect to cross-head displacement. 
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Table 4; For a given exit speed (V 2 = 5.23 mm/sec) and an initial thickness ti = 0.6 
mm, estimates of time spent under the roller bite during plastic straining. 


Plastic Strain 

h 

(mm) 

Time 

(seconds) 

0.05 

0.57 

0.33 

0.1 

0.54 

0.47 

0.15 

0.51 

0.57 

0.20 

0.48 

0.66 

0.25 

0.45 

0.74 


erature. The designed hxture Il60l . provides support to a long tail of the peel-specimen 
and eliminates any spurious effects due to gravity. When a stack of six layers is roll- 
bonded then a total of five bonded interfaces are formed. Peeling is done at the central 


interface. Figure 22 shows force versus displacement curve during the peel test. For 
all peel tests a cross-head speed of 15 mm/min was chosen. The steady-state peeling 
force P is used to estimate the rate of external work per unit advance of crack as 2P/b, 
where ‘b’ is the width of the specimen (typically 15 - 20 mm). In order to correctly 
determine the fracture toughness (Gc) of the plastically-welded interface, any amount 
of plastic work due to bending of peel arms must be subtracted from the total steady 
state work ll45l . Next, we present the mechanics of the peel test and methodology to 
estimate the correct interface toughness (Gc). The correction factor is adopted from 
||45l . The error bars in Gc (as shown in the main letter) are based on the variation when 
peeling force becomes steady. 


Elastica Analysis 


If the bending of the peel arm leads to only small elastic strains with possibly large 
rotations then it is referred to as an elastica. Figure [^schematically shows a peel test, 
where symbols have following meaning; 
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P is the vertical force applied by the upper-grip at ‘A’ [N] 

My is the moment applied by the upper-grip at ‘A’ on the peel-arm ‘OA’[N m] 

Mfo is the moment exerted at ‘O’ on the peel-arm ‘OA’ [N m] 
t is the thickness of the peel-arm [m] 

9 is the angle made by the tangent at any point along the peel-arm with respect to hor¬ 
izontal [rad] 

Oyj is the yield-strength of the material in tension [MPa] 
a,, c is the yield-strength of the material in compression [MPa] 

V is the Poisson’s ratio 
E is the elastic modulus [Pa] 

E’=is the plane-strain elastic modulus [Pa] 
t is the thickness of one peel-arm [m] 
b is the width of the peel-arm into the plane [m] 
s is the coordinate along the peel-arm [m] 

1= ^ is the moment of inertia of beam out of plane [m^] 
p = ^ is radius of curvature at any point on the beam [m] 

K = i is the curvature [m^ ' ] 

U is the elastic energy due to bending [J] 

Ge is the critical energy release-rate [J m^^] 
k= ^ is n constant defined for convenience [m^^] 


In Eigure 23 the crack tip location is marked at the location ‘O’. The upper peel- 
arm is shown as ‘O-A’ along with the forces acting on it. During the peel test, both the 
peel-arms are clamped in the tensile-tester. We shall carry out the analysis on the upper 
peel-arm (and symmetrical conditions hold for the lower peel-arm). 

When the width ‘b’ of the peel-specimen is considerably larger than the thickness; 
then plane-strain conditions prevail. In our case ‘b’ is 15 - 20 mm and ‘t’ is 0.2 - 0.4 
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Figure 23: Schematic of an elastica. For the sake of clarity only the forces and 
moment on the upper-peel arm at the crack tip due to lower arm are indicated as 
M^, and P. 


mm, hence plane strain bending scenario is assumed. In this plane strain problem is 
zero, and only and are non-zero stress components. Adapting elementary beam 
theory to plane strain, we have; 


^ 


(3) 


where, e^- = — ^ (y=0 is the middle thickness of the top elastica ‘OA’). The elastica 
is assumed to be inextensible, implying no stretching of the neutral axis under tension. 
The moment (M) at any section relates to the radius of curvature as: 
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( 4 ) 




E'l 

P 


Equilibrium of an element along the beam at point ‘B’, with coordinates (x,y) and 
angle 9, leads to 


dB 

E'I— = P{xmax-x) +Mf; 


(5) 


differentiating the equation |^with respect to ‘s’, give 


, d^Q dx 

E'l—^ = -P . 
ds^ ds 


(6) 


From the coordinate rule we know that. 


ds cos(0) = dx 


(7) 


and 


ds sin{Q) = dy. 

Substituting equation |7]in equation 


(8) 
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Pcos{Q), 


(9) 


E'l 


d^Q 

ds^ 


and on re-arranging, 


d^Q 

ds^ 


+ kcos{Q) = 0 


( 10 ) 


where, k—^. 


Now use a substitution ^ = v, and integrate the above equation to get 


— + ksm{Q) = Cl (11) 

where ci is the constant of integration, and it can be obtained through boundary 
conditions on v (^). 


Using equations 0 [n] we can find; 


S = 


X = 



dd 

-y/2(ci — ksm{Q)) 


cos{Q)dQ 
a/2(ci —ksin{Q)) 


( 12 ) 


(13) 
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Y = 



sin(9)dd 
■\/2(ci — ksin(Q)) 


(14) 


The above expressions also have an analytic solution in a specific case (as given in 
i45l ). It is important to note that Y and S can grow unbounded (depending upon the 
boundary conditions). 


We are interested in the conditions for: 


(i) Elastica, and 


(ii) Long peel-arm during steady-state 


The boundary conditions, in Figure 23 can be set as Mf = 0 at 0 = 7t/2, or curva¬ 


ture ^ = Kj, = 0 at 0 = n/2. 


Now using equation 11 we get 


Cl 



Thus, integrals [T^ [T^ [^can be reduced to 


Storat — 


rs rn/2 

ds = 

Jo Jo 


dQ 


\/2k{l — sin{Q)) 


(15) 
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^max — 


rx r 

/ dx = 

Jo Jo 


"/2 COi(0)t/0 

x/2k{l — sm{Q)) 


(16) 


Ymnx — 


ry f 

i- = / dy = 

Jo Jo 


V2 sin{d)dd 


y/2k{l — sin{Q)) 


(17) 


It is important to note that equations 15 and [^represent improper integrals of 
second kind since the function to be integrated is unbounded in the specified limits. 


This implies that elastica becomes vertical only asymptotically in the absence of 
any moment at the upper-grip. Numerical integration method can be employed by car¬ 
rying out the integration in the range [0, k/2). 


Next, we are interested in calculating the net energy stored in the (part-of) elastica, 
from 0 to 0, and this integral turns out to be a proper integral. 


For an element under bending as shown in Figure 24 we have 


e'i 

M= — = kE'I 
P 


(18) 


and 


a = 'Vids 


(19) 
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Thus, differential-energy stored in an element ds in bending from 0 up to an angle 


a (or curvature k) is 


^element — J element — J 


( 20 ) 



I I 
»ai 
I I 
I I 
I I 
I I 
I I 
I I 



Figure 24: Energy in bending. 


Using equations 


18 and 19 we can also write 


element (^) — ^ 


( 21 ) 


Thus, 


/■K 



k2 

E'Ik dK ds = E'I—ds 
2 


( 22 ) 
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Thus, the energy of the (part-of) elastica from ‘O’ up to any point along on it is 


given as 



(23) 


If converted in terms of 9 we get 



(24) 


For Qf — % 12, \]„et is the total energy of the elastica. The total bending energy 
is always finite even though the length of the elastica is unbounded, and it is worth 
emphasizing that this is correct when the elastica is assumed to be inextensible. We 
shall see some illustration graphs shortly. But, first we shall establish the conditions in 
which the limit of elastica is broken and plasticity starts due to bending. 

Conditions for Onset of plasticity 

If we assume isotropic material behavior, with Cy c = <3y t, the location of maximum 
curvature K occurs at ‘O’. Although glassy polymers may exhibit both kinematic and 
isotropic hardening, films with 42.3% PEG exhibited negligible changes in tensile yield 
stress even for plastic strain up to 25%; and therefore the original properties of films 
(before rolling) were chosen for analysis below. As an illustration. Figure shows 
the true stress-strain curves in tension for E15-alone-42.3% PEG films roll-bonded at 
different levels of nominal plastic strain. It is seen that yield points of films in tension 
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Figure 25: True stress-strain curves under tension for E15-alone-42.3% PEG, af¬ 
ter roll-bonded at different levels of nominal plastic strain. 

after roll-bonding at different levels of plastic strain are not much different than the 
yield point of the starting film. This indicates that effect of plastic strain during roll¬ 
bonding on the yield strengths of the films is negligible. 

The plane-strain condition and von-Mises yielding conditions imply that onset of 
yielding is marked by: 


^yield 

V 1 - V-l-V^ 

The maximum bending stress occurs at the base and is given as: 


(25) 


_ Mbt/2 _ 6Mb 
Omax - 


(26) 


Thus, setting <3max = we get the maximum bending moment at the base as: 

V l-V+v2 
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^b,max — 


^yield^^ 
6 \/l -V+' 


Since, Mj, — kE'I i.e. 


(27) 


^b^max — 


^yield^f 

6E'Iy/l — V+V^ 


(28) 


Now using equation El 


k2 

b,max 

2 


Pmax 


(29) 


This leads to: 


Pmax — 


<^yieldbt 


6E'{\ -v+v^) 


(30) 


If we stay within the limits of the onset of plasticity then elastica analysis applies. 
If the peel-arm is long enough then the steady-state bending energy of the elastica is 


constant (given by equation 241, and all the external work goes into the debonding 


process, and fracture toughness is given as: 


Gc = 


2P 

b 


The maximum Ge before the onset of plasticity is given as: 


(31) 


41 










Gr.max — 


^^yieldt 


6£'(1-v+v2) 


(32) 


Whenever experimentally measured steady-state peel force is greater than Pmax (es¬ 
timated by equation [30| ), a correction based on ll45l is applied. In steady-state the 
plastic work due to bending is subtracted from the total work and corrected interface 
toughness is estimated. 

As an illustration, let us consider b = 20 mm, t = 0.3 mm, E = 300 MPa, v = 
0.4, and Ge = nearly equal to 6 J/m^ (such that elastic limits are not crossed). The x- 
coordinate, y-coordinate and s based on numerical solutions are plotted in Figure 
Clearly, Xmax is finite whereas y-coordinate and total length, both, grow unbounded. 
However, the total bending energy of elastica up to some angle 0 are plotted in Figure 
27 As 0 approaches k/2—1 .57, the energy converges to a finite value. 


According to Figure we can say that the cut-off for the asymptotic behavior 
is at 60 mm i.e. 6 cm. This gives a sense that in a short length of 6 cm, steady-state 
elastica solution can be achieved. 


4.3 Lap-Shear Testing 

Preparation of lap specimens and shear-strength measurements were carried out in In- 
stron testing machine. A lap joint was assembled between two him layers, each layer 
being nearly 100 pm thick. The overlapping region was nearly A = 5x5 = 25 mm^ in 
area. A cross-head speed of 0.5 mm/min was chosen to apply desired compression load 
on the overlapping area. The sample was plastically bonded by pressing between two 
parallel (accuracy; 1 pm) Hats. Lap joint was tested for shear-strength in tension mode 
(at a cross-head speed of 15 mm/min). A snapshot of the test is shown in Figure 
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Figure 26: X,Y and S as function of angle 0. 



Figure 27: Bending energy U„e, as a function of angle 0. 



Figure 28: Lap shear-strength test specimen in tensile tester. 


43 













(■ -\ 



Displacement (mm) 

v_ ! _ ! _/ 


Figure 29: Force versus displacement curve during lap shear-strength testing. 
For this specimen a nominal plastic strain of 9.0% was imposed and the shear- 
strength (Oi=FmtH/A) was estimated to he 0.07 MPa. 


The peak force before failure divided by the bonded area was taken as the lap shear- 


strength. Figure 29 shows the force versus displacement during a lap shear-strength 


measurement. 


4.4 Mechanics of ‘hydrostatic die’ 


Piston 



Split-Bushing 

Cylinder 


I/— Specimen 


Plug 


33 


Figure 30: A set-up to achieve near hydrostatic compression. 

As discussed in the main letter and Supplementary video S2, the purpose of design¬ 
ing a ‘hydrostatic die’ was to explicitly show the role of active plastic deformation in 
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achieving sub-Tg, solid-state, plasticity-induced bonding. Figure 30 shows the CAD 


model of a ‘hydrostatic die’. Such a setup is capable of generating large levels of hy¬ 
drostatic pressure, while strongly limiting the plastic flow to negligible levels, when 
a circular stack of film with a radius equal to the internal radius of the cavity is com¬ 
pressed inside the die. 

We present a couple of analyses to demonstrate the principle of the ‘hydrostatic 
die’. Illustrations related to deformation theory of plasticity, as presented here, are 
borrowed in parts from lIlSl [^ . In what follows next, a boldface letter is to used to 
indicate a tensor variable. 


4.4.1 Elasticity Analysis 

First, we consider axisymmetric elastic compression of a film stack placed in the die. 
Here, all strains are assumed to be elastic and frictional forces are assumed to be absent. 
The solution to this problem is derived from the standard procedure of stresses in a 
thick-walled cylinder with a zero internal radius. A cylindrical coordinate system (r, 
0, z) is used. The principal stress components are denoted by Or, Ct© and a^, and the 
associated strains given as e^, e© and e^. All other shear components are zero. Due 
to axisymmetry and wall constraints inside the die we have Or = Ct© (= O, say) and 
Er = 6© = 0. Using the boundary constraints with the stress-strain relation of linear 
elasticity: 

1 -f V V 

E= -tr(a)I, (33) 

we And, 


(l-2v)(l+v) 


(34) 


Evez 

(l-2v)(l-fv) 


(35) 
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where, E and v are Young’s modulus and Poisson’s ratio, respectively. The stress 
tensor in terms of principal directions ir, e© and is given as o = o + O e© 0 e© 

+ Oj gj 0 e^. O can be decomposed into deviatoric part (o') and hydrostatic part (0^1, 
with Gm = {Or + O© + Oj)/3 denoting the mean normal stress), and written as o = o' + 
OmI. In the present scenario: 


3 


o = 


/ 20+g; 

3 


— 


V 


0 

0 


O-Oj 

3 


0 

0 

3 


0 

2o+o. 

3 

0 


2o+a- 

3 


The von Mises stress (Oy) is given as: 


Ov = 



(36) 


Thus, in the presence of die the von Mises stress (Gy die) and hydrostatic pressure 
(Pdie = -Om) are given as: 


Pdie — 


£|ez 


3(1-2v) 


(37) 


^v.die — 


1+V 


(38) 


In contrast if we imagined axisymmetric, unconstrained elastic compression, with¬ 
out any frictional effects then we would have o^ = and 0^ = 0© = 0. In such case 
the von Mises stress (Gy^„o-die) and hydrostatic pressure (pno-die) would be given as: 
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(39) 


Pno—die 


£|£zl 

3 


0—die — \^z 


(40) 


We emphasize that elasticity analysis is valid only up to the onset of plastic defor¬ 


mation. However, if we are within the elastic limit then equations 37 and H^show that 


large hydrostatic-pressure can build up during compression inside the die. Particularly 
in the limit as V —^ 0.5, pdie °°- 

In the video S2, a maximum compressive load of 40 kN is applied on a film-stack 
with radius 0.5" = 12.5 mm; corresponding to = —78.98 MPa. According to equa¬ 
tion |3^ if o, = —78.98 MPa, v = 0.4, and E=300 MPa then e, works out to be —0.12. 


Substituting je^| = 0.12, £ = 300 MPa and v = 0.4 in equation 38 Ov.rfie = 25.7 MPa. 


Clearly, Oy^die thus obtained is larger than the yield strength of the film. This indicates 
the possibility of plastic deformation even in the presence of the die. Next, we present 
an analysis based on the incremental (or “flow theory”) of plasticity, which accounts 
for plastic deformation. 


4.4.2 Incremental (“Flow Theory”) of Plasticity 

Here, we take into account the plastic deformation and demonstrate how little amount 
of plastic straining occurs when a film stack is compressed in the presence of the ‘hy¬ 
drostatic die’. 

According to the total deformation analysis, the total strain increment tensor (d^) 
is the sum of the elastic strain increment tensor (t/e'^) and the plastic strain increment 
tensor (deP): 


de = -I- deP 


(41) 
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The increment in elastic strain tensor can be derived using equation 33 and written 


— —d{tr{a))l (42) 

E E 

Under multi-axial loading the behavior of ductile materials can be described by the 
Levy-Mises equations, which relate the principal components of strain increments in 
plastic deformation to the principal applied stresses. In the present scenario this can be 
expressed as: 


P.P 




£ P 


(43) 


On the grounds of axisymmetric compression, similar to what was discussed in 
the previous section, we have = c® (= c, say), and therefore = e® (e, say). If 
we consider small intervals of time dt and call the resultant changes in normal plastic 
strains as r/ef = dt^ (= dt’’) and , it follows that: 


dt’’ = o'dl. 


(44) 


dtP = a[dX (45) 

here, dX is an instantaneous non-negative constant of proportionality which may 
vary throughout a straining programme. We further define following quantities to aid 
this illustration: 


|o'| = Vo': o' 


(46) 


a = 




(47) 
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= V d^P ; d^P 


(48) 


dzP = 


-dzP ; d^P = 



The flow rule can be written as; 


(49) 


dBP _ o' 
\d^P\ |o'| 


(50) 


If we assume that plastic flow is incompressible then the total increment in plastic 
strain can be written as; 


dz^ = —^dE^Cr^er— ^t/8^e0(8)e© + c/e^e2(8)ej (51) 

It is worth noting that at any instance neither the directions nor the relative mag¬ 
nitudes of plastic strain components change, and therefore we can write equation 
as; 




(52) 


By choosing D = — ®er— 0 e® + Cz 0 e^, we can re-write equation 


1 - 
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as; 


eP = 


(53) 


In equation 53 is negative during compression. Using the flow rule from equa¬ 


tion 


50 and the fact that increments in plastic strains are proportional to the principal 


directions (which are constant throughout the deformation history); we can rewrite flow 
rule in terms of total plastic strain at any instant as; 
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_ o' 

|eP| “ |o'| 

If we assume plastic deformation to continue at a constant yield Y, then; 


(54) 




(55) 


It is worth mentioning that the yield strength of polymers is usually a function of 
hydrostatic pressure and plastic strain (and the flow stress increases with increasing 
hydrostatic pressure and plastic strain, thus an assumption of constant yield stress is an 
under estimation). 

i we can write: 


By combining equations 53 54 and 55 


o = 


e? 2 


YD 


(56) 


Since ef is negative during compression, equation 57 can be re-written as 


a' = --YD (57) 

3 

The elastic strain increments, as given in equation occur both due to deviatoric 
stress and hydrostatic stress. At the onset of plastic flow (and continued plastic yielding 
at constant flow stress) the the deviatoric stress becomes constant (given by equation 
Hzl, after which there is no further contribution to elastic strains due to deviatoric stress 
components. However, the normal stress and the hydrostatic part of the elastic strains 
continue to increase. Thus, total elastic-strain can be written as: 




From equation 41 the total strain can be written as: 


(58) 
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(59) 


e = e^D-(i±^)(|^)D + 


Om A-2V 

3 E 


Now imposing the constraint that total strains in the r and 0 direction are zero i.e. 
Er = e© = 0 then equation [59| implies: 


Om E 4 7(1+v) 

3 1 - 2v ^ 2 3£ 

Thus, the overall stress tensor can be written as: 


(60) 


a 



E 

2(1-2v) 


e^I- 


y(i+v) . 

3(1-2v) 


( 61 ) 


If we choose values same as in the previous section, we get jef | nearly 2%. We em¬ 
phasize the fact that in polymers the plastic straining is accompanied with hardening, 
and the yield stress increases with mean normal pressure, therefore, je^ | = 2% is quite 
an overestimation. Lastly, the ratio ^ = 0.96 which suggests that the state of stress 
inside the die is ‘hydrostatic’. 

These calculations demonstrate that despite the large hydrostatic pressures there is 
only small plastic straining, due to which no bonding outcome is noted. 


4.5 Mechanics of Axisymmetric Upseting 

Unconstrained compression of film stack with initial thickness much smaller than the 
radius of the stack qualifies as a classic case of an upsetting problem. As shown in 
the Supplementary video S2, a film stack with an initial thickness of 0.84 mm requires 
peak loads up to 40 kN (which equates to a peak nominal stress of 78.9 MPa, much 
larger than yield strength of the polymer) in order to achieve a final thickness of 0.70 
mm due to compression. This can be attributed to frictional forces acting on the top 
and the bottom surfaces during compression. 

Figure [^schematically shows the upsetting of a film stack. Figureshows the 
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stress components along with the frictional forces acting on a cylindrical element. Now, 
we estimate the loads required to achieve plastic deformation in upsetting scenario and 
their comparison with the experimentally noted loads. 



Figure 31: Axisymmetric Upsetting of laminates. 


z 


t 




Figure 32; Stresses and frictional forces acting on an element during axisymmet¬ 
ric upsetting. 


Upper Bound Analysis for Load Estimation 


Upper bound analysis is one of the methods to estimate the deformation load and 
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the average forming or forging pressure. The analysis presented here is borrowed from 
1^ 13^ . The usual assumptions are: 

1. The deforming material is isotropic and incompressible. 

2. The elastic deformations of the material (and tool) are neglected. The tool is 
essentially rigid, and material behavior is perfect-rigid-plastic. 

3. The inertial forces are small and are neglected. 

4. The frictional shear stress, T, is constant at the die/material interface and is de¬ 
fined as X = fa = ma/y/3. 

5. The material flows according to the von Mises flow rule. 

6. The flow stress and the temperature are constant within the analyzed portion of 
the deforming material. 

In addition, following steps are key to invoking the upper bound analysis: 

7. Description of a family of admissible fields (using parameters to be determined 
later); these must satisfy the conditions of: incompressibility, continuity, and velocity 
boundaries. 

8. Calculation for the energy rates of deformation, internal shear, and friction shear. 

9. Calculation for the total energy rate and minimization with respect to unknown 
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parameters of velocity field formulation. 


One of the simplest solutions to upset forging assumes that a decrease in specimen 
height is compensated by an increase in width or radius without any bulging or barrel¬ 
ing. This implies that a rectangular solid deforms to a thinner and wider rectangular 
piece. In an axisymmetric scenario, cylindrical solids retain their cylindrical geome¬ 
try. Since the lines parallel to the height axis remain parallel, a parallel velocity field 
is assumed. A similar assumption on the velocity field is also made here, along with 
homogeneous deformation. 


Velocity and Strain Rates 

In order to estimate loads, first we need to calculate the velocity field and strain- 
rates. During upsetting, as shown in Figure we have already assumed that the 
volume is constant during the plastic flow i.e. the volume of the material moved in the 
z direction is equal to that moved in the radial direction (shortly we shall also impose 
plastic incompressibility during calculation of strain rates): 

nr^Vo = 2nrVrh 


or. 


Vr = Vor/^h (62) 

In the z-direction, can be considered to vary linearly while satisfying the bound¬ 
ary conditions at z = 0 and z = h. In the tangential direction, 0, there is no flow (due 
to symmetry). Thus: 
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( 63 ) 


Vz = -Voz/h 

From the Figure the increase in strain in 0 direction, i.e., the length of the arc, 
is given by: 


c/e© 


{r + dr)d®-rd@ _ dr 
rd& r 


(64) 


thus, the strain rate is 


c/e© dr 1 Vr Vo 
dt dt r r 2h 


The other strain rates are: 


(65) 


. _ ^ _ -Vp 
^ dz h 


( 66 ) 


tr = 


dVr 

dr 



(67) 


jrz = 


/8vv dvz \ 
\ dr J 


= 0 


( 68 ) 


j0Z = Y0r = 0 


(69) 


Thus, the effective strain rate is: 


e = 



(70) 


The strains can be obtained by integrating the strain rates with respect to time, i.e.: 
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E; 


(71) 


or with dh = —Vodt: 


Z 



Vodt 

h 


Ez 




Similarly, the other strains can be obtained as: 


E© — Er 


2 ho 


Ez 

2 


The effective strain is 


(72) 


(73) 


e = 


(74) 


Upper Bound Analysis 

In the upper bound analysis, the load is obtained by equating the rate of work done 
by the tool with the upper bound estimate of energy expended due to deforming mate¬ 
rial: 


If Et is the total energy rate expended due to material deformation, then Ej =LxV 
Qo), here L represents the forming load and V represents the velocity of the die. 

The total energy expended in material deformation itself can be expressed as sum 
of energy rates for deformation (Ed), internal shear (Es) and friction (Ef)' 


Et = Ed+Es + Ef 


or. 
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( 75 ) 


Et= aedV + / x|Av|^^s+ / XiVids 
Jv Jss JSF 

where, v is the relative velocity between the two zones of material when the velocity 
has internal shear surfaces and x = a/>/3; S indicates surface (internal or at die/material 
interface), v, is the die material interface velocity in the “i” portion of the deforming 
material and x, = m,d/v^ which is the interface shear stress at the “i” portion of the 
deforming material, m is chosen as 1 for the upper bound analysis here. 

The velocity field for homogeneous upsetting has been calculated before and a 
constant flow stress (c) is assumed. The deformation energy rate is given as: 

Ed = f otdV = (76) 

Jv h 

Es (internal shear energy rate) = 0, because there are no internal velocity disconti¬ 
nuities in the assumed homogeneous velocity held. 

The friction energy rate is given as: 


Ep =2 XiVids 
JsF 


where v, is the radial velocity given by equation 62 and ds = 2%rdr. Ep includes 
the friction energies on both the top and bottom surfaces of the deforming part. Thus, 


= fVjr 

Jo 2h 2h Jo 


(77) 


r, with X = mo j >73, 


p 2 a Vd . 
Ep = -nm—-—R 
3 V3 /z 


(78) 


Thus, total energy rate is given as: 
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Ej 


kR^oVd + -Tlm—-^R 
3 v3 « 


o Vd, 


3 



Figure 33: Strain-rate sensitivity of E3/E15 in 1:1-42.3% PEG. Nominal strain- 
rates are listed. 

and, the load is estimated as: 

Et 7 f 2 R\ 

+ (79) 

If we set R — 12.7 mm, h = 0.1mm, c = 10 MPa, we estimate the L = 40.414 kN. 
As shown in the upsetting video S2 (unconstrained case), tg = 0.84 mm, tf = 0.7 mm 
and therefore e = ln{to/tf) 0.18. Duration of the loading is approximately 30 seconds, 
and we use a flow stress value a = 10 MPa (for an estimated strain rate of 0.006 sec^^) 
from Figure |^(which shows strain-rate sensitivity measurements in tension). 
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